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SIMM  ARY 

' he  continuity  equation  has  been  applied  to  fatigue  crack  models  in  which  all  cracks 
begin  together  to  form  a  family  o  f  smooth  trajectories  proceeding  at  different  rates.  When 
the  rates  at  a  particular  crack  length  are  log-Normally  distributed,  it  is  possible  to  estimate 
the  effect  o  f  the  mean  and  variance  of  crack  life  in  a  reliability  situation.  For  a  logarithmic 
variance  0  1  of  crack  rates,  the  standard  deviation  of  life  is  reduced  by  14" „  approxi¬ 
mately. 

This  crack  model  has  been  applied  to  a  previous  one-crack  analysis  to  allow  random¬ 
ness  in  both  initiation  and  cracking. 
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Thi •  continidtx  equation  has  been  applied  to  a  fatigue  crack  model  in  which  all 
track  s  Begin  together  in  /arm  a  /ami lx  <>/  smooth  trajectories  proceeding  at  different 
rates.  When  the  rales  at  a  par  in  alar  crack  length  are  log-Sormally  distributed,  it  is 
possible  to  estimate  the  elicit  an  the  mean  and  utrutnee  of  crack  life  in  a  reliability 
situation  far  a  lagarithmit  uiriamc  0  1  ot  crack  rates,  the  standard  deviation  of 
life  is  reduced  hi  14",.  approximated 

This  model  has  hccn  applied  to  a  preuous  one-craik  analysis  to  allow  randomness 
in  bath  initiation  and  <  rat  k  mg 
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NOTATION 


Alt 

a,  (alt) 

00  =  0^0 
f  —  f(o\t) 

f 

g »  h 

G,  G(ta ) 

H,  H(ta) 

Ma 

Ma(u) 

Mt(u) 

M^(u) 

r(a)  —  r  (t  „) 
ra  =  r(0) 

R,  R(a,  t) 

Ro(a)  —  Rn(ta) 

Ra  =-  /?o(0) 

I 

ta 

t.\ 

la  =  I -la 

Ta 

It 

X 

x 

<ho 

<f>A(lA ) 

bra 
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•> 
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Random  crack  length  conditional  upon  time  t 
Crack  length  (at  time  t) 

Length  of  initial  cracks 
Crack  length  density 

Probability  density  function  according  to  subscript 
Heuristic  constants,  Section  4.3 
Defined  by  Equation  (4.5) 

—  T{1  a)  Ut„ 

Moment  generating  function  of  median  cracking  life 

Moment  generating  function  of  crack  life 

MGF  of  initiation  time 

MGF  of  total  life 

Risk  or  hazard  rate 

Hijack  risk 

Generic  cracking  rate 

Median  growth  rate 

Initial  median  growth  rate 

Overall  time  or  life 

Median  crack  growth  time 

ta  with  random  cracking 

Initiation 

Median  life  to  runaway  cracking 
Transform  variable  for  MGFs 
Random  factor  for  particular  cracks 
Particular  value  of  X 
Density  of  life 

Density  of  crack  time  with  random  rates 
rth  moment  of  median  crack  life 

—  In  .v 

Variance  of  .v 

Variance  of  median  crack  life 


1.  INTRODUCTION 

Until  now,  the  general  theory  of  structural  fatigue1,2  has  treated  the  growth  of  cracks  in 
structures  with  two  distinct  stages — viz.  crack  initiation  or  damage,  followed  by  crack  growth. 
The  idealisation  implied  in  those  papers  is  that  variability  is  introduced  by  random  initiation 
time,  while  structural  interactions  are  related  to  the  growth  of  cracks.  In  reality,  the  second 
assumption  is  fairly  accurate  but  there  is  randomness  in  the  crack  growth  exhibited  by  nomin¬ 
ally  identical  structures.  So  far,  allowing  for  this  has  defied  analysis  for  multiply-cracked 
structures;  the  present  report  extends  the  previous  work  to  allow  random  cracking  rates  in 
single-crack  models. 

It  begins  with  an  exact  treatment  of  the  crack  growth  models  of  Payne3  or  Hooke4  which 
is  then  combined  with  random  initial  lives.  Previous  results  for  life  and  moment  generating 
function  are  generalised  and  there  is  some  discussion  of  the  low  order  effects  of  random  crack 
rates. 

2.  DISTRIBUTION  OF  CRACKING  LIFE 

This  replaces  the  reliability  type  of  life  density1  previously  found  for  one-crack  models 
with  deterministic  cracking,  and  indeed  it  will  transpire  that  the  present  result  is  obtainable 
from  that  model  by  simple  substitution.  Before  beginning,  however,  the  randomness  in  crack 
rate  must  be  described  more  closely. 


2.1  Nature  of  Random  Cracking 

For  the  earlier  deterministic  cracks,1  the  possible  trajectories  were  a  group  of  curves 
obtainable  from  one  another  by  translation.  These  were  indexed  by  the  initiation  time  t0,  but 
all  satisfied  the  rate  da/dt  =  R(a),  a  being  crack  length  and  t  time,  while  R  is  a  single-valued 
function. 

In  practice  R  is  random,  as  in  the  first  trajectory  of  Figure  I .  Each  growth  increment  during 
cycling  is  random  and  there  is  also  an  overall  materials  effect  acting  for  every  load  applied  to  a 
given  structure.  The  first  of  these  two  types  of  randomness  corresponds  to  some  stochastic 
differential  or  difference  equation  and  leads  to  the  second  order  diffusion  term  in  the  Fokker- 
Planck  equation.  The  second  overall  type  of  variability  allows  smooth  crack  growth  for  any 
particular  structure  but  variations  between  them  as  shown  on  the  right  of  the  figure. 


a 


time 


FIG.  1  NATURE  OF  CRACK  GROWTH 


If  diffusion  is  ignored  but  overall  variations  are  allowed  for,  the  fatigue  model — though 
different — is  still  deterministic  and  follows  the  previous  continuity  equation 


Df 

Dt 


V 

ci 


R(a,  0  . 


cf 


ca 


=  — /(div  R(a,  t )  I  r(a )) 
for  crack  length  density  /(a  t)  given  time  t. 


(2.1) 


2.2  Rate  Divergence 

We  now  assume  that  for  a  crack  beginning  at  to,  dajdt  =  XRo(a)  which  replaced  R(a,  t ) 
above.  Here  X  determines  the  overall  crack  rate  and  is  a  function  of  a,  t  and  to-  Without  know¬ 
ledge  of  a  t,  X  is  random  with  density  fx(x),  log-Normal  where  necessary;  the  particular  trajectory 
corresponding  to  any  x  is  a  transformation  of  X  to  A  t. 

Then  the  divergence 
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with  /  and  /«  fixed. 


For  any  crack  daidi  integrates  to 
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(2.3) 


This  is  illustrated  by  Figure  2. 


2.3  Crack  Length  Density 

With  the  present  model,  the  characteristic  equations  of  (2.1)  are 

~df  =dt= 
f(xRaXa) + r(a)  -f- 1  /(/ — /o))  .v^o(a) 


(2.4) 


Along  the  crack  trajectories  (for  which  x  is  constant) 

-rflog  f^Ro'(a)da+  + 

R0(a)  xRoia)  x(t-t0)Ro(a) 


da 
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d  log  R0(a)  +  da  f  —  —  •  - . 

xRo(a)  taRo(fl) 


Integrating, 
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Now  let  ra  be  the  growth  time  of  the  median  crack  (jc=  I)  and  change  the  crack  length 
variable  a  to  ta.  the  time  for  the  median  crack  to  reach  the  same  size. 

Then 


log  f  =  A  —  log  Ro(a) 


_  I  f '• 

*Jo 


r(c)dc— log  ta- 


(2.5) 


This  is  singular  at  initiation  ta  =  0  because  all  cracks  here  begin  at  the  one  point  (/0,  oo) 
and  (2.5)  indicates  an  infinite  density  corresponding  to  the  concentrated  probability  that 
a\t0  =  an-  However,  the  equation  suggests  that  we  express /in  terms  of  a  modified  “density”  fta- 
Thus, 


/(«  t)f a 


A  I  \'a 

*J  o  ,lr 


r(r)dr)  . 


(2.6) 


where  A  is  constant  along  each  trajectory,  as  is  x.  Consider  another  median  growth  time  xt  for 
which 


fxr  —  „  -  exp  I  — 

Ro(xt)  1  1 


/  I  f  VT  J  \ 

I  -  I  r(c)dc  I . 

V  A  J  0  / 


Substituting  for  A  in  (2.6)  leads  to 


..  „  f(a  t,  x)xt  /  I 

f(a  t)t„  =  Rq(xt)'  exp  - 

R<Ata)  \  x 


r(c)dc 


(2.7) 


XT 


Here,  of  course,  /«  =  x(t~to)  and  since  x  is  constant  along  a  characteristic,  it  may  be  can¬ 
celled  from  each  side  to  give 


AT)  exp 

where  Ro  =  /?n(0),  the  initial  growth  rate. 


K 


V| 


r(c)dc  1  . 


(2.8) 


2.3.1  Boundary  Conditions 

These  will  be  expressed  in  terms  of  f(at)ta.  Being  supported  only  at  (to.  an),  they  must  be 
regarded  asymptotically.  So  far  we  have  not  used  the  assumed  density /.v(.v),  nor  any  assumptions 
about  asymptotic  behaviour. 


3 


As  in  Figure  3  we  assume  a  finite  initial  crack  rate  (unless  *  =  0)  already  implied  by  our 
use  of  R0  =  /?o(0).  Together  with  continuity,  this  presumes  asymptotically  linear  initial  growth 
consistent  with  the  two-stage  fatigue  assumption  of  a  fixed  finite  initial  crack  length.  In  welded 
structures,  where  immediate  initiation  is  more  plausible,  the  same  behaviour  agrees  with  the 
assumption  of  a  pre-existing  flaw,  restricted,  however,  to  a  constant  size. 


FIG.  3  APPROACH  TO  BOUNDARY  CONDITIONS 


As  r  -  0 

a  xR nr 

and  the  asymptotic  crack  length  density 

f(a  t )  =  /v(  V)//?or, 

so  that 

r /(a  t )  -  fxixVRo 

and  (2.8)  becomes 
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(/  t„)Rn(x{l 
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■  'Jo 
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r(c)dc 


(2.9) 


2.4  Density  of  Crack  Life 

At  a  given  time,  the  total  rate  of  attrition  is 


Ml  hi) 


0(1 


/v(v) 

U  I  u)  Rail  a) 


r(r„)exp 


( :jr 


r(c)dc  I  da 


(2.10) 


remembering  that  median  growth  for  time  t„  .v(r  in)  produces  the  crack  (o,  i).  Again,  ta 
is  the  natural  variable  of  integration  and  we  obtain 
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mh) 


f.xVaKt-to)) 

t-l0 


r(ta)  exp 


dta 


(2.11) 


When  cracks  are  small  one  expects  this  to  be  consistent  with  attrition  at  the  constant  risk  r o; 
we  show  that  this  is  r o. 

Near  t  =  to,  the  reliability  type  density  from  (2.1 1) 


r(ta)  exp  (. .  .)  ->  r0  exp  (-r0(/-t0)) 


(2.12) 


for  all  cracks  from  (to,  ao),  i.e.  all  values  of  x.  This  factor,  the  desired  asymptotic  behaviour, 
may  be  taken  outside  the  integration.  The  remaining  integral  is  asymptotically  one  since 


fx(tal(t-t0))  =  fx(x). 


a  density  function. 


3.  RUNAWAY  CRACKS 

These  will  be  treated  by  telescoping  as  in  Reference  1.  The  present  assumption  is  that  the 
median  crack  “runs  away”  when  ta  =  Ta  so  that  other  cracks  are  runaways  after  the  actual 
times  TalX.  Mathematically  this  means  that  r(Ta+ )  ->  oo  so  that  the  risk  integrated  past  Ta 
is  also  infinite. 


3.1  Reduction  of  Life  Range 

Equation  (2.1 1)  may  then  be  written 

‘  00 

Ta 


it 

1 

'Ta+ 

Jo  J 

Ta  —  * 

fx(x)  .  .  /I 
t  to  \  * 


rdc  I  dta 


and  we  first  consider  the  third  integral  in  the  limit  as  risk  uniformly  increases. 
Thus,  let  R  >  r(r„)  >  p  when  ta  >  Ta  so  that 

t'la 


t-t  o 
fa 


R.(t-to)  •'  "|  r(c)dc  >  p(t-t o)  (x=tal(t-to)) 

0 


in  which  case 

e-RM  t 0>  .  exponential  term  <  e~iAt  ta' 
Let  ./lv(.v)/(r  —  to)  be  bounded.  We  require  an  upper  bound  for 


fxre  h  "<• J 


*•  f"  /,( '■-) 

t  to  J  Ta-S-  '  ^o/ 


lo> 


dta 

t—to 


J  Ta  \ 

<  Re  r'T.  (.’{1  — Fx(Ta/(t —tg))}. 

If  now  there  are  finite  constants  C,  d  such  that  R  <  Cpd  then 

ftp  -  iMm -to)  .  /.!/■  /ul 

->0  as  p  ••>  oo 

Thus  if  Cprf  >  r(ta  ta  >  Ta)  >  p  the  third  integral  in  (3.1)  vanishes  as  r(Ta  +)  ->  oo. 

3.2  Telescoping 

Now  consider  the  second  term 
,  i'Ta  i 
I  —  ta 


H-J^e^ppTrjj 


r(c)dc  dta 


(3.1) 
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rdcj  dta  -  OUT),  (3.2) 

using  the  mean  value  theorem. 

Apart  from  the  factor  l/.v  the  integrand  above  is  a  probability  density  of  life  appropriate 
to  a  risk  r(c).  To  develop  this  idea  consider  a  particular  crack  trajectory,  along  which  .v  is 
constant.  Then  reliability  theory  indicates  the  conditional  life  distribution 

'iMt  /u,  -v)  —  1  —  exp  I —  I  r{x{t'  t{\)W 


Jx(Tgj(t  —  tp)) 
t  to 


'Ta 


r(ta)e\ p 


Ta 


1  exp 


r(  (■)(/<• 


0 


with  c  =  x(t'-to)  and  t„  —  x(t— r<>).  This  has  the  density  derivative 


(3.3) 


hit  to,  v)  =  xr(t„)  exp  -  r(c)dc 

x  I  0 


fa 


(3.4) 


the  integrand  of  (3.2)  apart  from  .v,  here  equal  to  (t-tu)ITa  ■  0(A T).  Thus  (3.3)  provides  the 
integral  in  (3.2)  and  with  runaway  cracks  (3.3)  is  one  if  ta  Ta.  With  all  substitutions  the  limit 
of  (3.1)  is 


4>{t  /o)  = 


fx(x) 

t-ta 


r(ta)  exp 


'j  “r(c)dc\dt, 


(3.5) 


In  this  generalisation  of  (2.11)  Ta,  the  median  time  to  runaway  cracking,  is  effectively  a 
constant  parameter.  The  second  term,  the  density  of  runaways,  is  similar  to  but  not  the  same  as 
the  corresponding  term  with  deterministic  cracking:  the  same  idea  of  concentrated  crack  life 
attrition  still  applies. 


4.  LIFE  DISTRIBUTION 

From  Reference  I,  the  life  distribution  for  translated  cracks 

[To 

4>(t)  r»e  r«'{l  |  <’  rnl«fit»)dF„(t  -  to) 

Jo 


/  fTn 

re  r"H  T-\1o(t  T„)  exp  I  —  I  r{c)dc 

where  if  r„  t  t„ 


dF„jdt„  -  r(t„)  exp 

from  reliability  theory  applied  to  crack  life. 


(4.1) 


(4.1  A) 


6 


The  moment  generating  function  here  is 


MAu)  -  M,(u  r0)M„(u)  (M,(u  r„)  I)  (4.2) 

u  —  r*u 

where  A//(t/)  and  Mtt(u)  generate  moments  of  initial  life  and  crack  time. 

It  Ma  is  replaced  by  another  function  reflecting  the  outcome  of  attrition  on  randomly 
cracking  structures  from  t  =  0,  the  unique  inversion  of  (4.2)  indicates  that  <j>(t)  is  similarly 
generalised  by  generalising  the  reliability  density  (4.1  A),  the  inverse  of  the  new  Ma(u). 

For  the  different  but  deterministic  crack  model  above,  this  replacement  density  is  (3.5), 
which  already  includes  any  runaway  cracks  so  that  the  third  term  in  (4.1)  is  subsumed  into  the 
second. 

For  completeness,  the  generalised  life  distribution  will  be  found  directly.  For  this,  one  needs 
J'(a't,  to)  which  is  (2.9)  in  the  present  instance.  The  formula  (3.5)  is  also  the  density  of  Ia  =  r— fo, 
the  crack  part  of  the  total  life.  Equation  (3.5)  then  indicates  that  rA  is  independent  of  initiation, 
as  one  would  expect,  and  the  total  rate  of  attrition  is  thus 
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where  from  (3.5) 
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(4.3) 


(4.4) 


Because  renewals  affect  only  the  initial  life,3  the  form  (4.2) 
when  there  are  inspections  or  renewals. 


may  be  immediately  applied 


4.1  Moment  Generating  Function  of  Crack  Time 


This  is  the  Laplace  transform  of  (4.4)  which  in  general  is  difficult  to  simplify.  To  investigate 
the  effect  of  small  departures  from  translatable  cracking,  we  shall  now  assume  that  X  is  log- 
Normal.  Then,  for  .v  near  I,  we  shall  try  to  relate  the  extended  MGF  Ma  to  that  of  the  median 
crack  Ma. 

From  the  Laplace  transform  of  (4.4) 
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A 


In  both  terms  put  ,v  c%  t a  -■  t«c  f  with  the  Jacobian  J  —  1.  In  the  first  term  ,  the  upper 

limit  of  .v  -  T„!t a  from  which  we  need  £.  On  the  right,  *e  have 

-v  T„i’itt„  leaving  ta  =  T„ 

with  this  upper  limit  for  any  value  of  f.  Thus,  the  range  of  $  is  ( -  oc,  x)  whilst  0  ta  <  Ta- 

In  the  second  term,  x  =  T„!ia  for  all  1. 1,  and  with  this  we  find  ta  --  tAx  =  Ta:  that  is,  the 

range  of  ta  is  just  the  point  T„.  On  the  other  hand,  this  identity  and  the  fact  that  0  tA  < 
means  that  .v  has  the  same  range.  Thus,  in  the  runaway  term,  the  two  substitutions  have  changed 
an  integration  over  1. 1  to  an  averaging  over  (.  It  would  be  possible  to  avoid  some  of  this  manip- 
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ulation  if  ihe  MGF  were  originally  defined  as  a  Stieltjes  integral  with  a  concentrated  probability 
of  runaway  cracking. 

With  these  changes 

p  „  j*x 

MAu)  rVak  '• 

J  0  J  x 


where 


(4.5) 


(></„)  m„ 


r(  <■)(/(■. 


4.2  Crack  Life  MGF  for  Small  Variance 


Now  assume  that  £  -  A(0.  a-)  and  expand  about  $  0.  In  this  section  and  the  next  steps 

to  successively  cruder  but  more  practical  approximations  are  marked  by  equations.  If  we  use 

G(i„)  ~  lit, j  ■  H(t„) 

and  take  three  terms  of  Gexp  f,  then 

|X 


M.M 


r(r„k 


Utii 


•'*  exp 


o 


<j2G-/2  \ 

,  '  }( I  -  <r2G)  1  Vfa 

I  CT-U 


.  c  nr,i  exp 


/  <r-6  2  \ 

(, .  <;  ,2- 


(4.6) 


It  is  now  convenient  to  subsume  these  two  terms  into  one  Stieltjes  integral  over  a  semi- 
infinite  range  with  the  probability  density  i IF„(I„).  The  terms  in  (4.6)  containing  variance  expand  as 


exp 


ij-G'~;2 
I  r  n-G 


I 

(I  o-G) 


I  •  H)  j<7'(]//2-  H3  ■  >,H')  ■  .  .  . 

ul„ |l  .(<r-(//-  3 H  I)  |<j'(  3 H  4\H -  3 H'  J H*)  ■ 

UFta- 1 1  IF-  3//  2)  j<r'<3  iH1  3 H:i  ■  \H^)  ■  .  .  .} 


0 (t/:,/„:'  6). 


(4.7) 


Here  the  sum  of  expectations  of  the  first  term  on  each  line  obviously  form  Ma(u),  generating 
moments  of  the  median  crack  life.  If*/  0.  the  kernel  corresponding  to  the  first  line  would  (if 
complete)  be  the  probability  density  (3.5)  of  our  present  model.  For  M , i,  the  expectation,  how¬ 
ever,  is  taken  over  median  crack  lives;  the  first  term  is  one.  so  that  the  variance  terms  on  the 
first  line  have  zero  expectations.  In  (4.7)  therefore,  multiples  of  the  variance  terms  of  the  first 
line  may  be  subtracted  from  the  others  within  an  accuracy  0((7(i).  This  simplifies  the  approxi¬ 
mation  to 


M.t(u)  ^  M„{u )  u 


I  2H )  i-r »(’(//-  //)  H-'))ilF„(i„) 


!, il¬ 


ia1  !"'-’( I  H)  "'(j1  H:,)\<IF„U„) 


(4.8) 
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where  the  MGF  of  the  median  life 

j*x 

Ma(u)  =  e  u'»dFa(ta). 

Jo 

Obviously,  the  t/,  u2  terms  in  (4.8)  represent  increases  in  moments  caused  by  the  Normal 
scatter  of  crack  rates.  Similarly,  the  first  line  of  (4.7)  can  be  interpreted  as  the  corresponding 
change  in  crack  life  density.  Although  all  the  integrals  converge  the  errors  introduced  in  (4.6) 
by  the  quadratic  approximation  for  exp— f  will  be  large  unless  o  <  1  which  is  true  for  fatigue. 


4.3  Heuristic  Approximations  of  Crack  Life  Moments 

In  this  section,  we  derive  a  rule  of  thumb  to  estimate  the  effect  of  crack  rate  variability. 
To  proceed  further  from  (4.7)  or  (4.8)  it  is  necessary  to  avoid  higher  moments  of  Fa.  These  are 
therefore  replaced  by  the  corresponding  cumulants. 

We  could  also  bound  the  risk  as 

H(ta)  g>ai/Jo  —  hi  a,  say,  (4.9) 

and  find  strict  bounds  for  the  moments.  This  amounts  to  approximating  the  given  crack-life 
risk  problem  by  a  constant  risk  problem.  Since  these  bounds  are  likely  to  be  very  conservative, 
we  shall  abandon  rigour  henceforth  and  regard  g  —  0(1)  as  an  averaging  correction,  retaining 
all  signs.  This  will  provide  heuristic  estimates  of  changes  and  indicate  their  possible  importance. 
Treating  (4.9)  as  a  simple  substitution 

"  M„{lt)-u\\a2fla  -  a2h{  i  —  jo2)j.l>a  —  ■  UV/MaJ 

-r  $W2K'(I  -  \o2)nia  —ho2fl3a  —h^Fia}- 

Now  substitute  cumulant  formulae  for  the  moments  //ra  with  a  view  to  then  neglecting  those 
above  k>  (we  retain  the  notation  //„,  aa2  for  the  first  two  cumulants).  This  is  done  because 
cumulants,  unlike  moments,  can  decrease  with  order.  Thus,  neglecting  higher  cumulants  and 
terms  in  a6  and  higher,  we  find 

Ma(u)  -  M„(u) 

-£<t2(I  -jo2(l  g))  •  IgV]  r  {oa2ln„)[ -  ga2(  1  - 1<72) - 2  ■  25gV  -  3gV] 

•  lg3oA°a(oalfltt)3} 

r  1  -lo*)-gc?-g*o*)-*aH2g^  9gV  +  1 5gV(a0//,B)2]}  (4.10) 

where  g  %  I .  If  we  now  make  all  the  gs  one,  this  eventually  leads  to 

ft  a! Fa  -  I  t0'6<7-(l  ;  cr2)— (<T„/^a)2[l  —  I  •5<ri(l  4  (<7a//*a)2)] 

(oaM2  -  I  —  er-( 5  —  CT- (2  -  25  :  {<r„in„)2{  14 -3(aa(/<(1)2}])  (4.11) 

where  a2  is  the  logarithmic  variance  of  crack  growth  rates,  while  //„  and  <t„  refer  1o  reliability 
type  crack  lives  based  on  the  median  crack  growth  curve. 


4.3.1  Example 

Equation  (4.11)  has  been  derived  from  reliability  based  lives  for  structures  which  crack 
immediately.  For  the  general  two-stage  model  of  fatigue  it  indicates  the  effect  of  random  rates 
on  cracking  times. 

The  model  of  crack  rates  is  the  coarse  randomness  in  which  the  possible  crack  trajectories 
are  randomly  scaled  versions  of  the  same  median  curve.  As  a  typical  example,  put 

a  0-2303,  onZ/ia  (c'r'i— I)  =  0-2589 

which  corresponds  to  a  log-Normal  life,  with  logarithmic  standard  deviation8  01  (to  base  10). 
To  base  10  the  crack  rate  standard  deviation7  is  also  0- 1 ;  the  quantities  in  (4.1 1)  refer  to  natural 
logarithms. 
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After  substitution,  we  find  ha;  Ha  1  -  0  03344,  (oAjaa)i1  ~  I  — 0 •  2561 7  and  oAjaa  = 
I  -0  1 3754. 

The  general  effect  indicated  seem  to  be  a  small  reduction  in  mean  crack  time  and  a  possibly 
more  important  reduction  of  its  variance. 


5.  CONCLUSIONS 

The  methods  of  Reference  I  have  been  applied  to  one-crack  models  in  which  all  the  crack 
trajectories  are  smooth,  but  form  a  family  in  which  the  growth  rate  for  a  given  crack  length  is  a 
random  factor  of  the  median  rate.  The  application  is  straightforward  and  it  is  simple  to  allow 
for  runaway  cracks.  There  is  no  simple  relation  to  the  results  for  median  crack  lives,  although 
asymptotic  effects  on  the  moments  are  suggested. 

Application  to  a  typical  example  seems  to  indicate  that  the  mean  life  with  random  cracking 
is  slightly  less  than  that  from  the  median  crack,  but  there  is  a  more  important  reduction  of 
variance.  These  results  are  being  sharpened  by  further  study  of  the  asymptotic  moment  generating 
function. 

The  model  described  above  refers  to  structures  in  which  cracking  begins  immediately.  In 
Section  4  the  generalised  density  corresponding  to  this  has  been  used  to  extend  the  previous1 
two-stage  model.  This  allows  randomness  in  both  the  initiation  and  crack  growth  stages. 
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